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Abstract 

We investigate energy (momentum) distributions in J ftp photoproduction and J/tp 
production in B meson decay. In particular the upper endpoint region of the spec- 
trum is examined where the effect of soft gluon emission from the cc pair becomes 
important. Constructing a model which is consistent with the so-called shape func- 
tion formalism we consider these fragmentation effects and show that the relevance 
of possible colour octet contributions in the photoproduction channel is still incon- 
clusive. 



*Talk given at the conference "Hyperons, Charm and Beauty Hadrons", Valencia, June 2000. To be 
published in Nucl. Phys. Suppl. 



1 Introduction 



Non-relativistic QCD (NRQCD) [|l|] provides the framework for a successful description of 
many heavy quarkonia production processes. Performing a systematic expansion in the 
typical velocity v of a (anti)quark inside the quarkonium NRQCD includes modes where 
the quark antiquark pair (QQ pair) is produced in a colour octet state and hadronize to a 
colourless quarkonium by soft gluon radiation. 

Like other effective theories NRQCD implies a strict factorization between the short- 
distance production of the QQ pair and its long-distance fragmentation into the quarko- 
nium. While one can compute the partonic process perturbatively to definite order in 
the strong coupling constant a s the non-perturbative physics has to be parameterized 
into so-called NRQCD matrix elements (0^( 2S+1 Lj)). They give the probability for the 
hadronization of a colour singlet/octet (c = 1/c = 8) QQ pair with spin S, orbital angular 
momentum L, and total angular momentum J into the quarkonium H and scale according 
to power counting rules of NRQCD, i.e. for example that the colour octet matrix elements 
are suppressed at least by factor of v 2 with respect to the leading order colour singlet 
contribution. 

Although NRQCD matrix elements describe the hadronization adequately in most 
quarkonium processes, the theoretical prediction of its energy distributions needs more 
effort: Towards the endpoint of the spectrum the quarkonium takes more and more en- 
ergy of the incoming particles, i.e. the phase space for radiating off soft gluons from the 
quark antiquark pair becomes less and less. Therefore one should expect that the colour 
octet contributions where the QQ pair must emit a soft gluon to get rid of its colour are 
suppressed for large values of the quarkonium energy. However, in the inelastic J/ip pho- 
toproduction channel naive NRQCD calculations yield a steep raise of the colour octet 
contributions in the upper endpoint region contradicting the observation of a rather flat 
spectrum ||. Thus we parameterize the fragmentation process by functions, the so-called 
shape functions, rather than by numbers to account for the phase space effect mentioned 
above. 

Technically speaking the necessity of shape functions instead of matrix elements is 
caused by the breakdown of the NRQCD velocity expansion near the endpoint of quarkonia 
energy spectra ||. When the quarkonium carries a fraction (1 — e) of its maximal energy, 
where e is a small number, the kinematic limitations from the hadronic process appear 
in the short-distance calculation as expansion in v 2 /e. While for small quarkonia energies 
(e v 2 ) the naive use of NRQCD is appropriate the endpoint region (e ~ v 2 ) needs a 
resummation of the leading twist terms 0((v/e) k ). Thus one gets a shape function for each 
production channel separately which smears the partonic spectrum. As result the endpoint 
moves from the partonic value given by the mass of the QQ pair to the hadronic one set 
by the quarkonium mass. 

The outline of this article is as follows. First we construct a model for the shape function 
which is consistent with the shape function formalism of [|J and models the region where 
e <C v 2 . Thereafter we apply our model to the decay B — > J/ipX and fit the model 
parameter A to the CLEO data. Assuming universality we finally transfer the results from 
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Figure 1: Diagrammatic representation of a general quarkonium production process. 
J /if) production in B decay to the photoproduction channel.fj] 

2 Shape function model 

Starting point for our analysis is figure [I]. Here quarkonium production is considered as a 
two step process: the short- distance production of the QQ pair described by the partonic 
cross section <5"Qo[ n ] and its subsequent fragmentation into the quarkonium H. Hence we 
get 

{2i,?2ii «^- H s ?/(S i * <wi " i(o ■ Fn[i) (i) 

where the shape function 

/dk 2 d 3 k 
27(27)^ ( 2 ^( P + l -PH-k)Mk;PH,P) (2) 

parameterizes the emission of a soft gluon field with total momentum k and invariant mass 
k 2 from a QQ pair of quantum state n = 2S+1 l}/j > and with off-shellness I = Pqq — P where 
P 2 = ^vp?q. is a radiation function that contains the non-perturbative physics of 

the hadronization process. Before defining this radiation function we elaborate the phase 
space exactly. 

To this end we restrict ourselves to the case of a single massless hard particle in the 
final state and refer to the quarkonium rest frame defined by pn = P = rather than to 
the center-of-mass frame (P- m = 0). After integrating over the hard momentum px and the 
soft momentum k we switch to light cone variables for I to reproduce the shape function 
formalism of Q. We integrate out l + = Iq + l z , l\ and the azimuthal angular and rewrite 

^For a more detailed presentation of the contents of this article please refer to S by Beneke et al. 
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the lo integration into a integration over the soft gluon energy ko- Then the final result 
reads 



a/3 (/3 2 +fc 2 )/(2/3) 
„2 



(27F) = ^J^ J dk f.H n (P m ,PJ, Px ) ■ ^ _ a) (3) 

n (a 2 +fc 2 )/(2a) 

where the integration bounds are functions of 

a = P m+ - M H , p = P in - - M H (4) 

which again are defined in the quarkonium rest frame. 

The hard subprocess in — > QQ[n] + X is given by the flux factor T and the azimuthal 
average H n = J d(f>/(2n)H n of the partonic amplitude squared. While this part can be 
calculated perturbatively we have to model the radiation function $ n . We make the ansatz 

$ n (k;p H ,P) = a n ■ \k\ b -exp(-k 2 JA 2 n ) ■ k 2 ex V (-k 2 /A 2 n ) (5) 

where the exponential cut-off reflects our expectation that the typical energy and invariant 
mass of the radiated system is of order A n ~ rngv 2 ~ several hundred MeV for n = ^q 8 ^, 
3 Pj 8 \ 3 S ( f'\ Note that there is no soft gluon emission necessary in the colour singlet case, 
i.e. $ n collapses to the delta function 5 4 (k). The normalization a n is fixed as described in 
p. We choose the other parameters as follows (c = 1.5) 0: 

6[^ 8) ] = 2 , 6[ 3 P (8) ] = bfsf] = , (6) 
A[^ 8) ] = A[ 3 P (8) ] = A , A[ 3 SJ 8) ] = cA . (7) 



3 Momentum spectrum in B —> J/ipX 

Let us now apply the formalism to the J ftp momentum spectrum in the semi-inclusive 
decay B — ► J/tpX. The leading partonic decay process b —> cc[n] + q (q = {d, s}) results 
in a Jfij) with fixed momentum but the hadronic decay spectrum is modified by both 
initial and final bound state effects. The Fermi motion of the b quark inside the B meson 
is accounted by the ACCMM model in a simple but satisfactory way. The fragmentation 
effects are considered by our shape function formalism. Colour octet contributions are non- 
negligible because their suppression due to a scaling factor off 4 ~ 1/15 from the NRQCD 
matrix elements is compensated by a factor C 2 ^/C 2 ^ ~ 15 from the Wilson coefficients in 
the partonic process. 

Thus we take the partonic amplitude squared H n (m b ,2m c ) for n = 1 Sq 8 \ 3 Pj 8 \ 3 S^\ 
and 3 jSf^ at tree level and perform the substitution 

2m c -> Mjk) = ^M 2 J/ip + 2Mj^k + fc 2 (8) 
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Figure 2: Sum of colour octet modes dBR[ 8 ]/dpR to the differential branching ratio of the 
decay B — > J/ipX compared to CLEO data 0. 



to allow the cc pair being off-shell before it emits soft gluons. Then the J/ip energy 
distribution in b quark decay is obtained by folding the partonic result with our shape 
function model. Finally we use this result as input for the ACCMM model to include 
initial bound state effects. It assumes an isotropic Fermi motion of the b quark inside the 
meson and a Gaussian momentum distribution with a width pp of several hundreds MeV. 
Furthermore the ACCMM model keeps the kinematics of the b quark decay 'in flight' exact 
treating the soft degrees of freedom inside the B meson as spectator of mass m sp . 

CLEO has measured the J/ip momentum spectrum on the T(4S) resonance ]7|. Hence 
we boost our result from the B meson rest frame to their laboratory system. We then 
assume that the colour singlet contribution is dual to the exclusive mo des B -> J/i)K^ 
!■ 

The result for the branching ratio with J/i/jK^ subtracted is shown in figure 0. It is 
clearly seen that the effect of cc fragmentation is necessary to reproduce the data for an 
ordinary choice for the ACCMM parameters {pp = 300 MeV, m sp = 150 MeV). In other 
words A = cannot describe the data without pushing pp up to values far above the ones 
successfully employed e.g. in semi-leptonic B decays. If we fit the shape function model 
parameter A we obtain 

A = 300;% MeV (9) 

which excellently fits to our expectation A ~ m c v 2 ~ Aqcd from the NRQCD velocity 
scaling rules. The uncertainties in (P) come from a different weighting between the colour 
octet contributions and from varying pp from 200 MeV to 500 MeV. 
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94/95 * 96/97 prel. 
■ 97 prel. 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 3: J/ip energy spectrum in photoproduction for different values of the transverse 
momentum cut (p™* = 1,2,5 GeV): colour singlet model (dotted), naive (dashed) and 
shape function modified (solid) NRQCD calculation versus HERA data for p T > 1 GeV 
and > 2 GeV respectively. 



4 Inelastic J/ip photoproduction 

In this section we discuss the energy spectrum in inelastic J/ip photoproduction. Again 
we take the squared amplitude H n from the partonic calculation and perform the sub- 
stitution @ 2m c -> M cg (fc).0 In this case it is more complicated to incorporate the shape 
function model (|J) since the incoming particles distinguish a symmetry axis in the partonic 
process g + 7 — > cc[n] + g. We will keep exact its perpendicular l± dependence even though 
it is of higher order in velocity scaling. Next we employ the parton density function to in- 
clude initial state effect and integrate over the proton momentum fraction of the gluon and 
the transverse momentum p? of the J/ip to get its energy spectrum da(jp — > J/ipX)/dz 
where z = (pjM-Pp)/(p 7 -p p ) is the photon energy fraction transferred to the J/ip. 

Finally a comment on the normalization of our result. The relative weights of the 
different production channels are set by the corresponding values of the NRQCD matrix 
elements f|. The overall normalization is fixed by adjusting the shape function modified 
curvature to the one from naive NRQCD in the region 0.1 < z < 0.4. Thus we account for 
the effect that partonically the c quark mass @ is higher than the value m c = 1.5 GeV 
used in the fits for the NRQCD matrix elements. 

Figure ^| displays the result for pt > 1 GeV and pt > 2 GeV plotted against data 
from HERA ||. Due to our normalization procedure the shape function improved curves 

* There is a misprint in the original paper || : In each mode of H n the factor gl should be substituted 
by gf. The numerics, however, remain unchanged. 
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coincide with the naive NRQCD calculation for low values of z. However, they display a 
strong enhancement in the region of intermediate z before they fall down to zero for z — > 1 
as anticipated from general arguments. 

The explanation for this (unphysical) strong enhancement lies in the structure of par- 
tonic colour octet amplitudes. Since H n is proportional to (1 — z cZ )~ 2 for n = and 
n = 3 Pj 8 ^ where z c5 is the corresponding quantity to z for the cc pair we get large contri- 
butions from z C c —> 1- Unfortunately soft gluon radiation in the final state decouples the 
energy (as well as the transverse momentum) of the J/ip and the cc pair, i.e. including our 
model we sample the partonic rate for z ct - > z close to one even though the cut on p? sets 
an upper bound on z. To fix this problem we repeated the analysis with a higher cut 
(cp. figure [3]). Now the hunch disappears because z max is small enough to keep z c5 away 
from one. 

In conclusion one should note that it is impossible to neglect or even rule out colour 
octet contributions in the photoproduction channel from the observation of a flat energy 
spectrum without increasing p^. 
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